ON A FINITE ELEMENT APPROXIMATION OF THE STOKES 
PROBLEM UNDER LEAK OR SLIP BOUNDARY CONDITIONS 

OF FRICTION TYPE 



TAKAHITO KASHIWABARA 

Abstract. A finite element approximation of the Stokes equations under a 
certain nonlinear boundary condition, namely, the slip or leak boundary con- 
dition of friction type, is considered. We propose an approximate problem 
formulated by a variational inequality, prove an existence and uniqueness re- 
sult, present an error estimate, and discuss a numerical realization using an 
iterative Uzawa-type method. Several numerical examples are provided to 
support our theoretical results. 



1. Introduction 

We consider the motion of an incompressible fluid in a bounded two-dimensional 
domain with some nonlinear boundary conditions, specified as the slip boundary 
condition of friction type (SBCF) or the leak boundary condition of friction type 
(LBCF). These boundary conditions were introduced by H. Fujita in [7], and sub- 
sequently, many studies have focused on the properties of the solution, for example, 
existence, uniqueness, regularity, and continuous dependence on data, for the Stokes 
and Navier-Stokes equations under such boundary conditions. Details can be re- 
ferred to in [7] itself or in [18], [M], [17], and [1], among others. Similar types of 
nonlinear boundary condition, such as subdifferential boundary condition or Tresca 
boundary condition, have been reported in |13j . [5], and [3], among others. 

The frictional boundary conditions under consideration have been successfully 
applied to some flow phenomena in environmental and medical problems such as 
oil flow over or beneath sand layers and blood flow in the thoracic aorta. Such 
applications have been discussed in [Tl|, [19], [21], and [20]. In these works, the 
finite difference method is used for discretization, and theoretical considerations 
such as convergence are not addressed. 

On the other hand, few studies have focused on the theoretical analysis of numer- 
ical methods for these boundary conditions, even if restricted to the Stokes problem. 
For example, Li and Li |15j proposed a finite element approximation combined with 
a penalty method for the Stokes equation with SBCF. They proved the optimal or- 
der error estimate; however, they did not focus on a numerical realization of their 
finite element approximation. 

The purpose of this work is to construct a comprehensive theory of the finite 
element method applied to flow problems with SBCF and LBCF, including all of 
the existence and uniqueness result, error analysis, and numerical implementation. 
In doing so, herein, we restrict our consideration to the stationary Stokes equation 
in a two-dimensional polygon. 

The remainder of this paper is organized as follows. In Section 2, we review 
the results for the continuous problems described in [7]. Weak formulations by an 
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elliptic variational inequality for SBCF and LBCF are also presented. In Section 
3, we prepare the finite element framework using the so-called P2/P1 element, and 
state several technical lemmas. 

Section 4 is devoted to the study of approximate problems for SBCF. We propose 
the discretized variational inequality problem, proving the existence and uniqueness 
of a solution. In the error analysis, we first derive a primitive result of the conver- 
gence rate 0{h™^'^^'^'^^^^) under the H^^'^-H'^ regularity assumption with < e < 2. 
Second, we show that it is improved to 0(ft,™'"^^'^^) under the additional hypothesis 
of good behavior of the sign of the tangential velocity component on the boundary 
where SBCF is imposed. A sufficient condition to obtain 0{h'^), which is of optimal 
order when e = 2, is also considered. Finally, we propose an iterative Uzawa-type 
algorithm to perform numerical computations, and prove that the iterative solution 
indeed converges to the desired approximate solution. 

Section 5 is devoted to the study of approximate problems for LBCF, in a manner 
similar to Section 4. However, it should be noted that unlike in the case of SBCF, 
we have to explicitly deal with an additive constant for the pressure. As a result, 
sometimes there exist multiple solutions for the pressure, especially its additive 
constant; other times it is uniquely determined. Moreover, in an error analysis, we 
can only obtain the convergence rate 0{h™"^^'^/^-^/'^^), because of the error of the 
additive constant of the pressure. If we can eliminate the influence of this error, 
the same rate-of-convergence as in the case of SBCF is realized. 

In Section 6, several numerical examples are provided to support our theory. We 
observe that the results of our computation capture the features of SBCF and LBCF 
and that the numerically calculated errors decrease at 0{h?) for both. Section 7 
presents the conclusions and discusses some future works. 

The author learned about Ayadi et al. [2] after the completion of the present 
study. They treat the finite element approximation for the Stokes equations with 
SBCF, using the PI bubble/Pl element. Some numerical examples are presented, 
and an error estimate is announced without a proof. 

2. Settings and results of continuous problems 

2.1. Basic notation. Let be a polygonal domain in R^. Throughout this paper, 
we are concerned with the Stokes equations written in a familiar form 

-i^Au + Vp = f in divu = in Q, (2.1) 

where v > is the viscosity constant; u, the velocity field; p, the pressure; and /, 
the external force. As for the boundary, we assume that F := dQ is a union of two 
non-overlapping parts, that is, 

r = rouri, ronri = 0, 

where FcFi are relatively nonempty open subsets of F. Moreover, Fi is assumed 
to coincide with whole one side of the polygon for the sake of simplicity. Two 
endpoints of the line segment Fi are respectively denoted by Mi and M^+i; the 
meaning of these subscripts is clarified in Section [3?T| 

We impose the adhesive boundary condition on Fg, namely, 

u = on Fo, (2.2) 

whereas on Fi, we impose one (and only one) of the following boundary conditions 
of friction type: 

U„ = 0, \<Jt\<9, arUr + g\ur\ = 0, (2.3) 
called the slip boundary condition of friction type (SBCF), and 

Ur^O, |cr„|<g, (TnUn + g\Un\ = 0, (2.4) 
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called the leak boundary condition of friction type (LBCF). The function g, called 
the modulus of friction, is assumed to be continuous on Fi and strictly positive on 

ri. 

Here, the definitions of the symbols appearing above are as follows: 

n = *(?ii,rt2) = outer unit normal to the boundary Fi, 

T = ^{n2, —ni) = unit tangential vector to the boundary Fi, 

Un — u ■ n = normal component of u on Fi , 

Ur ^ u ■ T = tangential component of u on Fi, 

, , If dui duj \ p ^ . , 

Cijiu) = - — h = component of rate-oi-stram tensor (1 < 7 < 2), 

2 \0Xj oxi J 

Tij{u,p) — —pSij +2i>eij{u) = component of Cauchy stress tensor (l<i,j<2), 



2 

a(u,p) — Tjj {u, p)nj ^ — stress vector defined on Fi, 

cr?i = <7n{u,p) = normal component of stress vector defined on Fi, 
a-r — ariu) — tangential component of stress vector defined on Fi. 



Remark 2.1. (i) n and r arc constant vectors because Fi is a segment. 

(ii) (Tr does not depend on p, which is verified by a simple calculation. 

(iii) In (2.3) and (2.4), g acts as the threshold of the tangential and normal stress 
beyond which non-trivial slip and leak on Fi may occur, respectively. This is why 



the boundary conditions (2.3) and (2.4) are said to be "frictional.' 



2.2. Function spaces. We use the usual Lebesgue spaces L'^ (il) , L°° (il) and 
Sobolev spaces H^{n) for a nonnegative integer r, together with their standard 
norms and semi-norms (for a space of vector-valued functions, we write L^(ri)^, 
and so on). i7°(r2) is understood as L^(f2), and i?o(ri) denotes the closure of 
C^in) in We put Q = L^ifl) and 



Q = Ll{n) = |v e L^in) J ijjdx = 

H^(n) is also defined for s e R+ \ N by the norm 



E 

|q|— r 



f2xSl 



u _ ni\2+2e 



dxdy 



where a G is a multi-index and s = r + 6, where ?■ e N, < 6* < 1. 

We also use the Sobolev space H'^iT) defined on the boundary F for s > 0. 
i/°(F) is understood as i^(F), and we put 



L2(Fi)= 77eL2(Fi) 



rj ds = 



where ds denotes the surface element of F. The usual trace operator defined on 
onto i?^~^/^(F) is denoted hy tp i-^ V-'lr for s > 1/2; however, we simply 
write tp instead of V'lr when there is no ambiguity. Since n and r are constant 
vectors, we immediately obtain the following: 

Lemma 2.1. Let s > 0. For every (p E H^{r)'^ satisfying cp — Q on Fq and (p^ = 
{resp. (pr — 0) on Fi, it holds that 

HWH^ir)-^ = Il0r||ff=(ri) {resp. ||(/)||^.(r)2 = \\(pn\\H=iri))- 
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In addition, we require the so-called Lions-Magene space HQg{Ti) (see [121 Sec- 
tion 1. 11]) with its norm defined by 

Ni . =(\\vr^ +( \^ds^ 

where p{x) — dist(x, {Mi, Mm+i}) is the distance from a; e Fi to the extreme 
points of Fi along Fi. We use this space for only one purpose described in the 
following lemma. 

Lemma 2.2. The trace operator maps {%jj£H^{Vt) | ^^ = on Fq} onto _ffQg(Fi). 
Proof. See [101 Theorem 1.5.2.3]. □ 

Remark 2.2. The lemma implies that the extension to F by zero of an arbitrary 

1 1 
function in Hqq{Ti) belongs to H^{T). 

Now we let V = H^{n)'^ and introduce the following two closed subspaces of V: 

Vn^{ve H^{nf I w = on Fo, t;„ = on Fi}, (2.5) 

K = {w e iyi(f7)2 |i, = on Fo, t;^=OonFi}, (2.6) 

which corresponds to the velocity space for SBCF and LBCF, respectively. Com- 
bining the above two lemmas with the usual trace theorem, we see that 

Lemma 2.3. (i) For every v £ Vn [resp. v G Vt), it holds that 

^At, e Hoo(ri) and ll^rll^i^p^^ < C||i;||/i-i(n)2 

1 

{resp. i;„|ri e ffoo(ri) and ll"nll^i(r^) ^ ChWminy), 
with the constant C independent of v. 

(ii) Every rj £ Hqq{Ti) admits an extension v £ Vn {resp. v £ Vt) such that 

= 77 {resp. Vn = 77) on Fi and \\v\\m(^ny < C||?7||^i 
with the constant C independent ofrj. 

2.3. Bilinear forms and barrier terms of friction. Let us introduce 

a(w, w) = 2z/ V / eij{u)eij{v)dx {u,v £ H'^ {^f), (2.7) 

b{v,q)^- I divv qdx {v £ H\n)^ , q £ {Q)) , (2.8) 

J{V)= I 9\v\ds {j^£L\T,)). (2.9) 

The bilinear forms a and b are continuous with their operator norms ||a|| and 
respectively, being bounded. As a readily obtainable consequence of Korn's 
inequality ([TH Lemma 6.2]), there exists a constant a > such that 

(Vw e w = on Fo). (2.10) 

This implies that a is coercive on Vn and V^. We simply write j'("r) and j(w„) 
to express j{vT\ri) and j(fn|ri), respectively. Then, j{vr) and j{vn), called the 
barrier terms of friction, are continuous functional on V because j{vi) is bounded 
on L2(ri). 
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2.4. Green's formula. For all {u,p) G H^+'-^Q)'^ x iJ'(f^) with e > 1/2 satisfying 
divu = 0, we obtain Green's formula as follows: 

{—vAu + Vp,v)L2(^^^2=a{u,v) + b{v,p)— / a{u,p)-vds (Vw e 



where the stress vector a{u,p) is defined in Section 2.1 In fact, the line integral 
over r appearing in the right-hand side is well defined because cr{u,p) e iJ'^ 2 (r) C 
L^(r). However, if we have only a lower regularity, say {u,p) G iJ^(ri)^ x L'^{il), 
then the definition of a(u,p) in Section 2.1 becomes ambiguous. We thus propose 
a redefinition of (j{u,p) as a functional on iJj (r) as follows. 

Definition 2.1. Let (u,p) G H^{nf x L'^{Vl) with divu 0. When -vAu + Vp 
is represented by / G ^^(il)^ in the distribution sense, that is, 

we define ct(u,p) G {H2{T f)' by 



{v€H\nf). (2.11) 



Here and hereafter, for a Banach space X, we denote the dual space of X by X' 
and the duality pairing between X and X' by (•, 

Remark 2.3. The functional a{u,p) is well defined according to the trace theorem 
and the fact that the right-hand side of (2.111 vanishes if w = on F, i.e., v G 
iJg(r2)^. In addition, this definition of a(u,p) agrees with the previous one if u and 
p are sufficiently smooth to belong to 7J^+'^(ri)^ x H'^{Vt) with e > 1/2. 

In particular, we see that Ut and cr„ are characterized in H^^iTi)' by 

{(Jr.^Vr) 1 = a(u, u) -t- - (/,w)L2(f2)2 



a(M, - (/, w)l2(o)2 

respectively, in view of Lemma |2.3[ 



G Vn), 

[v G \^.), 



2.5. Variational formulation to the Stokes problem with SBCF. Herein we 
assume / G L'^{fl)'^ and g G C^(Fi) with g > on Fi. With Vn defined by (2.5) and 

o 

Q — -L§(f2), we introduce a weak formulation of the Stokes equations (2.1 1 under 
( p^ and ([53| as follows. 

o 

Problem 2.1 (PDE). Find {u,p) £ Vn x Q such that ar — crr{u) is well defined 
and the slip boundary condition of friction type (2.3 1 is satisfied, that is, 

(V« G Vn), 

(Vq e Q), 



' a{u,v) + b{v,p) - (crr,WT)L2(ri) = {f,v)mn)^ 
biu,q) = 

CTr/g G L°°(Fi) and |crr| < g a.e. on Fi, 

arUr + g\Ur\ = SL.C. ouFi. 



(2.12) 
(2.13) 
(2.14) 
(2.15) 



Note that G L^(Fi) follows from (2.141, and thus (2.12) makes sense. 
Another formulation by a variational inequality proposed in is 

o 

Problem 2.2 (VI). Find {u,p) G K x Q such that 

j a{u,v - u) + b{v - u,p) + j{vr) - j{ur) > if,v - u)L2i^^y {\/veVn), (2.16) 

\b{u,q) = (VgGQ). (2.17) 

The following theorem concerning the existence and uniqueness is essentially 
derived from Theorems 2.1-2.3]. 
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Theorem 2.1. (i) Problems PDE and VI are equivalent in the sense that {u,p) € 

o 

VnX Q solves Problem PDE if and only if it solves Problem VI. 
(ii) Problem VI has a unique solution. 

Remark 2.4. In [7 , another definition of = J't^ is employed and it is supposed 
that r is smooth, with Fq n Fi ^ 0. However, some shght modification, which is 
not essential, makes the proofs in [7| applicable to our own situation. 

2.6. Variational formulation to the Stokes problem with LBCF. As in the 



previous subsection, using Vr defined by (2.6) and Q ~ L (Q), we introduce a weak 



formulation of the Stokes equations (2.1) under (2.2) and i2AK as follows 



Problem 2.3 (PDE). Find {u,p) & V^- x Q such that cr„ = (t„(m,p) is well defined 



and the leak boundary condition of friction type (2.4) is satisfied, that is, 

a{u,v) + b{v,p) - {an,Vn)L^ri) = {f,v)L^(ny2 (Vu £ K), (2.18) 

b{u,q)=0 (VgeQ), (2.19) 

cr„/g e L°°(Fi) and |cr„| < a.e. on Fi, (2.20) 

CnUn + g\un\ = a.c. ou Fi. (2.21) 
Note that (T„ € L^{Ti) follows from (2.20), and thus (2.18) makes sense. 



Another formulation by a variational inequality proposed in 7J is 
Problem 2.4 (VI). Find {u,p) G K x Q such that 

a{u,v-u) + b{v-u,p)+j{vn)-j{un)>{f,v-u)L2(^a-^2 (VveK), (2.22) 

b{u,q) = (Vgeg). (2.23) 

We recall the existence and (non)uniqueness theorem derived from |71 Theorems 
3.1-3.3 and Remark 3.2]. 

Theorem 2.2. (i) Problems PDE and VI are equivalent in the sense that {u,p) € 
Vr X Q solves Problem PDE if and only if it solves Problem VI. 

(ii) Problem VI has at least one solution, the velocity part of which is unique. 

(iii) // {u,p) and {u,p*) are two solutions of Problem VI (therefore, Problem 
PDE), there exists a unique constant (5 € R such that 

p = p* + S and an{u,p) = <7n{u,p*) — S. 

(iv) Under the assumptions in (iii), if we suppose u„ ^ on Fi, then 6 = 0. 
Namely, a solution of Problem VI is unique. 

Remark 2.5. Although the definition of cr„ = —p+i/'^i^ and the hypotheses on the 
boundary in [7] are apparently different from ours, we can complete the proof with 
only a non-essential modification of the original one in [7] . 

3. Finite element approximation 

3.1. Triangulation. Let {.%i]h be a sequence of triangulations of a polygon 57, 
where h denotes the length of the greatest side. As usual, we assume that 

• n Tj is a side, a node, or for all T^, Tj ^ ST^ii ^ j). 

• I^Tj = 51 and the boundary vertices belong to F. 

i 

• When h tends to 0, each triangle in contains a circle of radius Kh 
and it is contained in a circle of radius K'h for some constants K,K' > 
independent of h. 

• Each triangle has at least one vertex that is not on F. 
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The one-dimensional meshes of T and Fi inherited from the triangulation are 
denoted respectively by S'h and ^hIti- °^ nodes, we use 

= set of all vertices of triangles in 

= set of all midpoints of sides of triangles in 

^o,h = To n 

ri,fc = Ti n = {Ml, M3/2, M2, • • • , M„, M„+i/2, M„+i}, 

f = Ti n = ri,ft \ {Ml, M„+i}, 

where the subscripts of Mj's are numbered such that 

• Mi^s, for i = 1, 2, • • • , m + 1, arc all vertices of triangles in which are 
located in Fi and are arranged in ascending order along Fi. 

• Mi_|_i/2 is the midpoint of Mj and Mj+i for i = 1, 2, • • • , m. 

In particular, To.h H Ti^h = Fq n Fi = {Mi,Mm+i}. Wc denote each side with 
endpoints Mi,Mj_|_i by Cj = [Mj,Mj_|_i] and its length by |ej| = |MiMj_|_i|, for 
I = 1, 2, • • • ,m. 

3.2. Approximate function spaces. We employ the so-called P2/P1 element, 
defining VhCV = H'^{Q)^ and C Q ^ L^{Vt) by 

Vh = [vh e C\nf I vhW e ^2(r)2 (vt e .%)^ , 
Qh = {qh e c°(n) I qh\T e ^i(r)' (vr e ^h)} , 

where ^k{T) denotes the set of all polynomial functions of degree on T (A; = 1, 2). 

o 

To approximate Vn, Vr, and Q = LKfl), we set 

Vnh = {vh e Vh I Vh{M) = (VM e To,h), Vhn{M) = (VM e Ti,h)} , 

Vrh = {vh e Vh I VhiM) = (VM e Tom), Vhr{M) = (VM e Fi,/,)} , 

Qh = Qh^Ll{n), 

together with 

Vh^Vhr\Hl{nf, 

o 

Vnh,a = [vh G Vh I b{vh,qh) = (Vgh G Qh)], 
Vrh.a = {vh G Vh I b{vh,qh) =0{yqh G Qh)}- 
Here, ii/j„ and Vhr denote Vh ■ n and Vh ■ t, respectively. By a simple observation 
we see that Vnh C Vn, Vrh C Vr, Qh C Q, Qh C Q, and Vh = Vnh n H^ifl)"^ = 

Vrhr\Hl(9)\ 

The quadratic Lagrange interpolation operator Xh : C°(fl)^ — ^ Vh and LP'- 
projection operator lih '■ Q ^ Qh are defined in the usual sense, that is, 

IhveVh and {Ihv){M) ^ v{M) (Vi; G VM G S,,), 
^hq&Qh and {q-Ilhq)qhdx = {^q € Q, \/qh € Qh)- 



It is easy to verify that IhV G V^/i (resp. IhV G if w G fl C"^(0)^ (resp. 

00 

V e Vr n C°(n)^) and that Uhq G Q/, if qh G Q. The following results for the 
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interpolation error are standard (for example, see 0) and are used without special 
emphasis in our error analysis: 

\\v-Ihv\\H^^n)-<Ch'-\\v\\Hi^.(a)- {yveH'+'inf), (3.1) 
\\q-'nhq\\mn)<Ch'\\v\\H^(n) {VqeH^m, (3.2) 

where < e < 2 and the constant C > depends only on il. Note that _ff^+'^(51)^ C 
C°(r2)^ for e > 0. Furthermore, the estimate on the boundary, together with Lemma 
|2.1| and the trace theorem, gives 



\\Vr - {Ihv)r\\L^^^) < C/l 3 1 1 1 1 ^ l ^ < C/l^ +1 W 1 1 l + . (o)2 

(vesp. \\vn - (Ihv)n\\L^r,) < Ch^+^WvnW^i^,^^^^ < Ch^+^\\v\\Hi+>(n)2 



(3.3) 
(3.4) 



for all veVnD H^+'{nf (resp. w G K n H^+\VLf). 

For approximate functions defined on the boundary Fi, we define 



A, 



[^Xh e C"(Fi) I ^ih\e e ^2(e) (Ve G <?ft|rj, ^^/.(Mi) =/ift(M„+i) =o}. 



|/ift(M)|<l(VAfefi,h)} 



A/i = e Ah 

By a simple calculation, we find that (see also Lemma [3.3[ i)) 

Aft = |wftr|ri w/i e = ji'/mlri e Vrhj ■ 



The space Aft also becomes a Hilbert space if we define its inner product by 

^ m 



i=l 



(Aft,/ift e Aft), 



(3.5) 



which approximates J^^gfih^h ds by Simpson's formula. Here and in what follows, 
we occasionally write g^. Aft i , • • • instead of ^(Afj), Aft(A/j^i ), • • • , and so on. 

o 

Since g is assumed to be positive on Fi (particularly, on Fi,ft), (■, Oa^ is indeed 
positive definite. Let us denote the projection operator from the Hilbert space Aft 
onto its closed convex subset Aft by Proj;^^ . It is explicitly expressed as 



-1 

ProjA,(M/.)(A'/)= <(/i/.(M) 
for each fih £ Aft 



if tih{M)>l 
if \^^h{M)\ < 1 
if fihiM)<-l 



(VAf e Fi,ft), (3.6) 



Finally, to approximate j given in (2.9), we introduce jh as 

jhiVh) = Q^\ei\(^g^\Vha\ +45,+ i|ryft,,+ i| + gi+i\vh,t+i\^ (Vh € Aft), (3.7) 

i=l 

again with Simpson's formula in mind. Clearly, jh is a positive, continuous, and 
positively homogeneous functional defined on Aft. This definition of jh is motivated 
by Section IV.2.6] and [8, Section n.5.4]. 

3.3. Inf-sup conditions. Hereafter, we denote various constants independent of h 
by C and those depending on h by C{h), unless otherwise stated. In this subsection, 
two types of inf-sup conditions concerning the approximate spaces of the velocity 
and pressure are considered. The first one is the "TJo-Lq" type and well known, 
while the second one is the "H^-L^^' type and seems to be new. 
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Lemma 3.1. (i) There exists a constant /3 > independent of h such that 

aw II ^ b{vh,qh) ^ ° 

PhhWmn) < sup — (Vq/i e Qh). 

(ii) Let fi and f2 be functions in i^(r2)^ and L'^{fl), respectively. Then there 

o o 

exists a unique {uh,Ph) ^Vh x Qh such that 

o 

a{uhiVh) + b{'"h,Ph) = (/i,w/i)l2(o)2 (Vw,, G Vh), 

O 

b{uh,qh) ^ if2,qh)L^{n) (yqh&Qh)- 
Moreover, {uh,Ph) satisfies 

+ \\Ph\\L2(Q.) < C(||/i|lL2(f^)2 + ||/2||l2(0)), 

where the constant C depends only on \\a\\,a,l3. 

Proof See [3 Chapter 12]. □ 

o 

Remark 3.1. Since C Ki/u we immediately deduce from (i) that 

All II ^ b{vh,qh) ,^ ^ ° . 

p\\qh\\L-^(n) < sup T — 7 lyqh e Qh)- 

Lemma 3.2. There exists a constant /3 > independent of h such that 

/3|k/.||r^(o) < sup ..^^^"'^^^^ [yqh^Qh). 

Proof. Let us take an arbitrary ph £ Qh and define r^^ G by 

1 



m — i + 1 



/ ■ _ m+2 m+3 2m+2 \ 

to/2 ^ 2 ' 2 ' ■ ■ • ' 2 J 



where ^ := l)L2(n) (iTil denotes the length of Fi). According to Lemma 

|3.3[ i), which is preceded by this lemma only for the sake of convenience, we can 
choose Uh € Vrh such that Uhn = rjh on Fi and 

Il"/i||_f/i(f2)2 < C'll'7/i||_H-i/2(ri)- (3.8) 
Then, by direct computation we deduce that 

Uhnds=[ ??/ids = -^^^ -(ph, l)i2(Q), (3.9) 

Fi JFi ^ 



and that 

WVhfH^r,) - + ^) < C\{pH, l)L2(n)P < C||p,,||i.(f,) 

The latter estimate implies 

WVhW^kir,) ^ C'h/.lkMF,) < C|K|U2(a). (3.10) 
From (3.8) and (3.10), we have 

\\uh\\HHnr < ^^11^/111^2(0). (3.11) 

For iih constructed above, it follows from Lemma 3.1 ii) that there exists a unique 

o o 

(K^Ph) e X Qh such that 

a{ulvh)+b{vh,pl) = (V«„eT4), (3.12) 

O 

b{ul,qh) = [ph, qh)L^n) - b{uh,qh) i^qh G Qh), (3.13) 
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together with the estmiate 

WKWn^ny < C\\ph + diviihWmn) < CiWphWmn) + \\uh\\HHny) 

<C\\ph\\mn)- (3.14) 



Here we have used (3.111 to derive (3.14). 

Now, setting Uh — Uh + u^^ E Vrh and decomposing ph as ph = pi + Sh, where 

o 

Ph '■= Ph - (Ph, i-)L^ (n)^/p \ e Qh and 5h := {ph,l)L^{n)^ /l^ {\^\ denotes the area 
of fi), we see from (3.13) and (3.9) that 

biuh,Ph) = b{'U'h,pl) + b{uh,Sh) ^ {Ph,pl)L^a) - / divu,i da; 



= iPh,pl)mn) - 5h / Uhnds = (ph.pDmn) + {Ph,S)L^n) 



= \\Ph\\L^n}- 

On the other hand, it fohows from ( 3.11[ ) and (3.14) that 

||wh||_f/i(o)2 < C\\ph\\L^{a). 
From (3.15) and (3.16), we conclude 

b{vh,p,-, 



(3.15) 
(3.16) 



sup 

This completes the proof. 



b{uh,Ph) II 



□ 



Remark 3.2. We can regard this result as a discrete analogue of [IF, Lemma 2.2]. 

3.4. Discrete extension theorems. Let us investigate some discrete extensions 
of functions given on the boundary Fi to that defined on the whole domain fl. 

Lemma 3.3. (i) Every rjh e admits an extension Uh G Vnh {resp. £ Vrh) 
such that 

UhT = Vh {resp. Uhn = Vh) on Fi and \\uh\\H^(ny < ■ (3.17) 

(ii) For all i% € A/j {resp. r]h £ A^Ci -Lo(ri))) ^e can choose Uh in (i) such that 

Uh e Vnh,a {resp. Uh e Vrh,a)- 

Proof, (i) Let rjh G A/j. We discuss only the construction of Uh G Ki/i, because we 
can construct Uh £ Vrh in a similar manner by replacing n with r and vice versa. 
Define a piecewise quadratic polynomial iph on F by 

>/.|ee^2(e)2 (Vee^^), 

MM) = (VA/ e Fo,0, 

>,.„(M) = and (l>hr{M) = Vh{M) (VA/ e Ti,h). 



We find that 

(j)h = on Fo, 
and thus we obtain 



Hn on Fi, 



and 



^r = Vh on Fi, 



(3.18) 
(3.19) 



with the aid of Lemma [2?T] 

Now according to the property of the discrete lifting operator (see [U Theorem 
5.1]), there exists Uh £ Vh satisfying 

= on F and \\uh\\HHn)^ < C\\(j)h\\ (3.20) 



(r)2 



We conclude Uh £ Vnh and (3.17) from (3.18)- (3.20). 
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(ii) First, take an arbitrary rjh e A/j and consider an extension to Vnh,(T- It 
follows from (i) that there exists Uh G Vnh such that Uhr ~ Vh on Fi and 

(3.21) 



For such Ufi, by Lemma 



3.1 



we can find {u'^,p 
a{ul,Vh) + b{vh,pl) = 

together with the estimate 



^Vh X Qh satisfying 

{yvh e 14), 

o 

{yqh &Qh), 



(3.22) 
(3.23) 



(3.24) 

where the last inequality holds f rom (| 3.21[). Now, choosing Uh ~ u'^ + Uh G Vnh, 

o 

r/ii because €V/i, and 



we deduce that £ Vnh,<j from (3.23), that Uhr 



UhT 



that \\uh\\m{ny < C||?7,,||^i ^^^^^ from (3.21) and (3.24|. 

Next, we let ?7;i G n -Lq(Fi) and construct Uh G Vrh in the same manner as 
above by replacing n with r and vice versa. Then, it remains only to show that 

o 

b{uh, 1) = because we already know that b{uh, qn) = if qn G Qh- We can verify 
b{uh, 1) as 

b{uh, 1) = - / divu,j dx = - I u^n ds = - rj^ ds = 0. 
Jn JFi JFi 

This completes the proof. □ 

3.5. Properties of (•, •)Ah and jh- Let us establish several relationships between 
the inner product of and the functional jh, given by (3.5 1 and (3.7), respectively. 
We use a signature function sgn(a;) in the usual sense defined by 

1 (a; > 0) 



sgn(a;) 







-1 



{x = 0) 
(x < 0). 



Lemma 3.4. (i) If Uh G Vnh {resp. Uh G Vrh) and Xh G A^, then 

{uhT,Xh)Ah < jhiuhr) {resp. {uhn,Xh)Ah < jh{uh„))- 

(ii) Under the assumptions o/(i), the following properties are equivalent: 

(a) {Uhr,Xh)Ah >ih{uhT) {resp. {Uhn,Xh)Ai, >3h{uhn))- 

{h) {uhT,Xh)Ah. = jh{uhT) {resp. {Uhn,Xh)Ah ^ .ih{uhn))- 

(c) {uhT,tJ-h - A/0 Ah < {resp. {uhn, l^h - Xh)Ah < 0) {Vfih G A^^). 

(d) If M GF i,ft and Uhr{M) ^ {resp. Uhn{M) ^ 0), then 

\h{M) = sgn{uhr{M)) {resp. \h{M) = sgn(u^„(M))). 

(e) Xh = Proj^JXh + puhr) {resp. Xh = Proj^^ (A^ + /OUh„) (Vp > 0). 

(iii) When Xh G Ah, the following properties are equivalent: 

(a) XheAh. 

(b) {r]h, Xh)A^< jh{Vh) {^Vh eAh). 

(iv) When Xh £ Ah, the following properties are equivalent: 

{a) {vh,Xh)A,^o {ym& ^hnLl{Ti)). 

(b) There exists a unique constant Sh G R such that 

Sh 



Xh{M) 



g{M) 



(VM GFi,/,). 
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Proof. We establish statements (i) and (ii) only for the case Uh G Ki/i, because the 
proof remains valid when Uh € Vth, with n replaced by r and vice versa. 

(i) This is obvious because \\h{M)\ < 1 for all M € Ti^ if Xh € A/j. 

(ii) (a)=>(b) Since we have already proved the converse inequality in (i), state- 
ment (b) immediately follows from (a). 

(b) =>(c) Let (b) be valid. Prom (i), it holds that 

{UhT, IJ-h - Xh)Ah = {UhT, IJ-h)hh - jhiUhr) < (V/ifc e Ah). 

o 

(c) ^(d) Assume that (c) is valid and consider an arbitrary M e Ti,h such that 
Uhr{M) ^ 0. Let us define fXh € by 



\Xh{N) if N€ri,h\{M} 

lsgnK,(M)) if N = M. 



When M e S^, we can write M = Mi for some 1 < i < m-|- 1. Now, by assumption 
we have 

iuhr,fJ.h - Xh)Au = ^^(|e,_i| + |e,|) {\uhr{M)\ - Xh{M)uhAM)^ < 0. 

This implies that A/j(M) = sgn(Mft^(M)) because \Xh{M)\ < 1 and uut^M) ^ 0. 
Similarly, when M e S^, we can write M = M^^i for some 1 < i < m. Then, by 
assumption we obtain 



{Uhr,f^h - A/,)a, = ^g{M)\ei\(\uhr{M)\ - Xh{M)uhr{M) 



<0, 



from which Xh{M) = sgniuhriM)) follows. 
(d)=J>(a) If (d) is true, then we see that 

{UhT,Xh)Au = Q^\ei\(^gi\UhT,i\ + ^9i+^\UhT,i+l \ + 9i+l\UhT,i+l\^ = jhiuhr)- 

i=l 

(c)<;4>(e) This is a direct consequence of a general property of projection opera- 
tors. In fact, we obtain 

(W'h,r. - Xh)K^ < (Vnh e A,,) 

<^ {Xh+ puhr - Xh,Hh - Xh)Ah < (V/^?i G Ah^p > 0) 

<^ A,, = Proj;^^ {Xh + puhr) (Vp > 0). 

(iii) (a)=>(b) This is already shown in (i). 

o 

(b)=>(a) Let (b) be valid and consider an arbitrary M e T\^h- Define r]h € Ah 



by 

rih{N) 



if N&Ti,h\{M} 

-hi or - 1 if N = M. 



When M G EJj, we can write M = Mi for some 1 < i < m + 1. By assumption, we 
obtain {rjh, Xh)AH < JhiVh), which leads to 



^(|ei_i| + \ei\^g{M)(^Xh{M) ± l) < 0. 



This implies that \Xh{M)\ < 1. We obtain the same result when M G in a 
similar way. Therefore, we conclude that Xh & Ah. 
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(iv) (b)=>(a) Let (b) be valid and consider such 6h- Because Simpson's formula 
is exact for quadratic polynomials, for all rjh G H Lq{Ti), we have 



i=l 

™ I e ■ I ^ 

4 = 1 

= Sh / rihds = 0. 



(a)=>(b) Let (a) be vahd and consider i G {2, 3, • • • , m}. Let us make r]h G A/; 
vanish except on ei„i U e^. Then, statement (a) is equivalently written as 
(a') For all 't]h.i-^iVh.ijVh,i+^ G if we assume 

|ez-i|(4%^4-i +Vh,^) + \ei\ir]h,r+4:Vh,i+0 = 0, 

then we have 

\e^-l\\^4,r|l^^,_l{gXh)i-^ +VhA9>^h)^] + \ez\[rihA9>'h)z + 'ir]i^,i+^{g>^h).,+ ^^ = 0. 

Now, if we take r]h\ei_iuei such that rjh/i. = and 4??,^ ^_ i |ei_i| = —^r]h,i+^ kil = 
1, it follows from (a') that {gXh)i_i = (.gA/i)j_|_ i . Similarly, if we take r]h\ei-iUei 
such that i_ 1 — and ?7/i_i(|ei_i| + \ei\) = — 477^^ i |ej| = 1, it follows again from 
(a') that {gXh)i = {gXh)^+i. Hence, {gXh)^-^ = {gXh)i = {gXh)i+^- 

Repeating the above procedure for i — 2,3, ■ ■ ■ ,m, we conclude that there exists 
(5h G R such that 

g{M)Xh{M)^Sh (VMGfi,/,). 
This completes the proof. □ 

The following mesh-dependent inf-sup condition is important to deduce the 
unique existence of the Lagrange multiplier Xh G A/^, which appears in Sections 
I3and[5l 

Lemma 3.5. There exists a positive constant Ph depending on h such that 

MriH\W< sup (Vr;„GA,0 (3.25) 

resp. phhhWK, < sup j,''^'";'^^^^'' (V7?;,GA,)y (3.26) 
vhi^v^h \\vh\\m(ny'- J 

Proof. Because both || • || 1 and || • Ha^ are norms defined on A/,, which is of a 

(Fi) 

finite dimension, they are equivalent. Hence there exists a constant C{h) depending 
on h such that 

hh\\^^^^^^<C{h:)hH\W (Vr/„G A,). 
Now, we let -qh & and choose Uh G Vnh satisfying ( 3.17[ ). Then we conclude 

sup 7^ — 1^ > 7^ — 1^ = — 1^ >C[li)- — ll^/hllA^ 

vh^v^h \\'>^h\\H^n)2 \\Uh\\H^n)2 \\uh\\H^ny^ 

>cmm\W- 



We can obtain (3.261 in a similar way. □ 
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3.6. Error between j and jh- We begin with some generalization of [9i Lemma 
IV. 1.3] concerning the error between j and jh, which is necessary later in our error 
analysis. 



Lemma 3.6. (i) There hold 

jhiVh) < C|l77/jL2(ri) (yVh G A,,), 
WVhh^ < C\\Tjh\\L^(ri) C^Vh e Ah), 

with the constant C depending only on g and Ti . 
(ii) If g ^ C^{Ti), then for all < s < 1, we have 

\jh{Tlh)-j{Vh)\ < Ch'\\T]h\\H^{r^) i^Vh e A,,), 
with the constant C depending only on g and Ti. 



(3.27) 
(3.28) 



(3.29) 



Proof, (i) Let r]h £ Ah- On each segment = [Mi, Mi^i], take two points denoted 
by -^i+i E^nd M^^^s, whose meaning is understood naturally, for i = 1,2, •• • ,m. 
Let us define a piecewise constant function rh{gr]h) on Fi by 

rhigVh) 

in 

• K. ^ ^ "^6 "'"e "'"6 -* 

(3.30) 



i=l 



where xa denotes the characteristic function of A C Fi. 
Then we have 

3h{rih) = I \rh{9Vh)\ds<\V,\^/^\\rh{gilh)\\ 



L^(T^)■ 



(3.31) 



By direct computation, it follows that 

i=l 



i=l 

m I I 

< \9? E ^ + K^h + '"v'+i) 

i—l 
i—l 

- Vh.z 'nh,i+i + 2'q^ .j^^i{r]ha + Vh,i+i)^ 

^l^^Plgl^WrihWh^r^)- (3-32) 
Here we have used the inequality 

+ 4y2 + < 2{x^ + z^) + 8y^ - xz + 2y{x + z) (3.33) 



4 y 



X + z 



+ -(x-z)2>0 (Va;,y,zeR) 



to derive the third line. We conclude (3.27) from (3.31) and (3.32 1. 
The estimate (3.281 follows similarly if we remark that 
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(ii) Let rifi £ A/^. First, from the proof of (i), we see that 
\jhiVh) - j{Vh)\ < / \rh{grih) ~ 9Vh\ds 



< |rir/'||r^(m)-mllL^(r,) (3.34) 

< \Ti\'/m\rh{grih)\\mr,) + hVhWmr^)) 

< CIMlhf,)- (3.35) 

Before giving an estimate of \\rhigrih) - grihllL^iVi) which involves |l?7/i|lHi(ri), it 
should be noted that if e ^2 (R) we have 

, / a + 



(x) — 4>{a) — [x ~ a) 



so that 















J a 





dx 



{a+b)/2 



\t-a\'\<j)'it)\Ut 

b 

< 2{b-af I \(P\t)\^dt. 



(3.36) 



In view of the Taylor expansion of g, we apply (3.361 to deduce 

fM. , 1 rM. , 1 



Il{e^) 



\gzVh,i-gvh\ ds 



Mi 



IgiVKz - giVh + giVh - gr]h\ ds 



Mi 



< 2 



•5 „2 



Mi 



9i iviut -Vh) ds + 2 {g, - g) 7?^ ds 



Mi.i 



Mi 



< 2sup|.gp • 2 



2 pM 1 

" I,,' |2 



< max 



sup|5p suplg'l 



|r/;pds + 2(^^sup|g'| 

\'^i\'^\\Vh\\HinM,,M.^i]) 



9 ' 18 

for i = 1, 2, • • • , m. By a similar discussion, we have 



/2(e,) 



M. 



M. 



-gVh?ds < Ch^\\r]h\\%n^^), 



M 



\gi+\Vh-i+\ -gvhl^ds < c/i^ II HI (e.), 



+4 



M. 



\g^+i7jh^i+i - gijhl'^ds < Cft.2||f?h||Hi(e.), 



2 r^i^i 

" \Vh? ds 

Mi 



(3.37) 

(3.38) 
(3.39) 
(3.40) 



for each i. Therefore, it follows from (3.34) and ( 3.37 1-( 3.40) that 

\ih{->lh) -j{Vh)\'^ < C\\rh{gi]h) -gVh\\l2(^r,) 



1=1 

m 



i=l 
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SO that 



|jft('7/0 - jirih)\ < C/i||?7/i||^i(ri) 



(3.41) 



As a consequence of (3.35) and (3.411, we obtain the desired inequahty (3.291 by 
Hilbertian interpolation (see Chapter 14]) between L'^{Ti) and H^(ri). □ 

As will be shown in Theorems |4.2| and |5.3| below, the leading term of the error is 



that between jh and j, which is estimated by (3.291 with s — 1/2. However, under 
some additional conditions, we can obtain a sharper estimate than (3.291. 

Definition 3.1. An element rjh e Ah is said to have a constant sign on every side 
if, for any i = 1, 2, • • • , m, either of the following conditions is satisfied: 
(a)?7h|e, >0 or (b) %|e, < 0. 

Remark 3.3. Let rj^ G A/j have a constant sign on every side. If r//j > on e^-i and 
?7;i < on Ci for some i, then iih{Mi) — 0. 

Lemma 3.7. Let g G C^(ri). If r/h G A/i has a constant sign on every side, then 

\jhiVh) - j{Vh)\ < Ch^Wrjhh^ri)- 

Moreover, if g is a polynomial of degree < 1, then jhivh) is exact, that is, 

jhiVh) ^ jiVh)- (3.42) 

Proof. Let rjh G A/j have a constant sign on every side. Because rjh > or rjh < 
on Ci for each z = 1, 2, • • • , m and g is positive on Fi, we have 



rhigVh)\ ds 
g\Vh\ds - 



rhigVh) ds 
gVh ds 



where rh{grjh) is defined as (3.301. Summing up these terms, we obtain 



„ fit „ IIL n 

jhiVh) = / \rh{gVh)\ ds^"^ \rh{gVh)\ ds = '^ I 

„ m „ m „ 

jiVh)= / g\Vh\ds = Y^ / g\r]h\ds = Y^ / grj^ds 
Consequently, it follows that 

m „ 

lihiVh) - jiVh)\ / (rhigVh) - gVh)ds 

,■1 J ei 



rh{gr]h) ds 



(3.43) 



Let gh denote the linear Lagrange interpolation of g using the nodes in S'^^ n Fi^/j. 
Namely, gh is continuous on Fi and affine on each side — [Mi, Mi^i], satisfying 
gh{Mi) = g{Mi) for i = 1, 2, • • • , m. Then the Taylor expansion of g implies 

\gh{x) - g{x)\ < y sup \g"\ (Vx G Fi). (3.44) 

Now, let us estimate each term appearing in the summation on the right-hand 
side of (|3l3| by 



(rfiigVh) - 9hVh)ds 



+ / \gh- g\\vh\ds. 



(3.45) 
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Since Simpson's formula is exact for cubic polynomials, we can express 
f 1^1™ 



rh[gVh)ds+ —(gh,^+^ - g,+ i)r]h^,+ i. 



Thus, due to (3.44), the first term of (3.45) is bounded from above by 

1 



12 



|e,|/i sup|5 I |?7v+i| 



Note that there holds (of. ( |3.33[ )) 

\Vh,^+^\^\(i^\ 



< 



,15 |e 



4 15 



for i = 1, 2, • • • , m. Then, the sum of the first term of (|3.45|) is estimated as 

(rhigm) - ghVh)ds 



E 



1=1 

<^h'snp\g"\ (^h,,+ if\e.\^ 
<l/i2sup|g"|-2h.|U.(r,)|ri|V2 



El 



(3.46) 



Next, the second term of (3.45) is estimated by sup \g"\ \rih\ ds, which gives 
E^ \gh^g\\Vh\ds<^hHvip\g"\J^ \rj,,\ds<^h^supW'\\\m\\LHr,)\ri\'/^ 



(3.47) 



Hence we conclude from (3.43), (3.46), and (3.47) that 

\jhiVh) - j{Vh)\ < Ch'^\\r]h\\L^^^) 



If g is a polynomial of degree < 1, then both terms of (3.45) vanish because 
gh — g, from which (3.42) follows. This completes the proof. □ 

4. Discretization of the Stokes problem with SBCF 

4.1. Existence and uniqueness results. We propose approximate problems for 
Problem VI (therefore. Problem PDE) in the case of SBCF as follows. 

o 

Problem 4.1 (VI^). Find {ufi,ph) € Vnh x Qh such that 

(a{uh,Vh-Uh)+b{vh-Uh,Ph)+jh{vhr)~jhiuhT)>if,Vh-Uh)L2^i-iy {yvhGVnh), (4.1) 

\b{uh,qh)=0 (Vg^eQh). (4.2) 

Problem 4.2 (VI^.o-). Find Uh E Vnh.a such that 

a{uh,Vh - Uh) + jhivh) - ih{uh) > {f,Vh - Uh)L^(n)^ i\fvh G Vnh.a)- (4.3) 



18 



T. KASHIWABARA 



Problem 4.3 (VE,i). Find {uh,Ph,^h) € Vnh x Qh xA/^ such that 

' a{uh,Vh) + b{vh,Ph) + (w/ir, A/i)a,, = ifiVh)L2(ny (Vv/i G K/i), (4.4) 

< b{uh,qh)--0 (yqheQh), (4.5) 

, {uhr,f^h - Xh)A^ < {Vfih e A,,). (4.6) 

Recall that we are assuming / G L^(ri)^ and 5 E C'^{Ti). We first establish the 
existence and uniqueness of these approximate problems. 



Theorem 4.1. (i) Problem VI^.o- admits a unique solution Uh G Vnh,, 
more, it satisfies the following equation: 



Further- 
(4.7) 



(ii) Problems Ylh,a-, VI;i, and VE/i are equivalent in the following sense. 

(a) // Ufi G Vnh.rr is a solution of Problem VI^^o-, then there exists a unique 

o 

Ph G Qh suck that {ufi,pf,) solves Problem VI/j. 

o 

(b) // (uh,Ph) G Vnh X Qh is a solution of Problem Ylh, then there exists a 
unique Xh G Ah such that [ufnPh, ^h) solves Problem VE/j. 

o 

(c) If {uh,Ph, ^h) G Vnh X Qh X Ah is a solution of Problem YE h, thenuh solves 
Problem Ylh.a- 

Proof, (i) Since the bilinear form a is coercive on Vnh and the functional jh '■ 
Vnh ^ R is convex, proper, and lower semi-continuous (actually, continuous) with 
respect to the weak topology, we can apply to Problem Ylnh,a a classical existence 
and uniqueness theorem for second-order elliptic variational inequalities (see [H 
Theorem 1.4.1]). Thus, there exists a unique Uh G Ki/i,cr such that (4.3 1 holds. The 
equation (4.71 follows from (4.3) with Vh = and 2uh. 

(ii) (a) Let Uh G Vnh. a be a solution of Problem Ylh.a- Taking Uh ±Vh as a test 



function in (4.3), with an arbitrary Vh G Vh H Vnh,ai we obtain 

o 

a{uh,Vh) = {f,Vh)mny (V-y/i G V), n Vnh,^)- 



Moreover, from Lemma 



3.1 



i), we deduce the unique existence of ph G Qh such that 

o 

a{uh,Vh) + b{vh,Ph) = (/, w/i)l2(s;2)2 (Vv/i G Vh) (4.8) 

by a standard argument. 

Now we let Vh G Vnh be arbitrary. It follows from Lemma 3.3 (ii) that there 
exists some Wh G Vnh,a- such that Wh — Vh on F, which implies 



and 



jhivhr) = .jhiwhr)- 



Vh-WhE Vh 

Since Uh,Wh G Vnh.a, we conclude from (|4.3|), (|4.8|, and (|4.9|) that 



(4.9) 



a{uh,Vh - Uh) + b{vh - Uh,Ph) + jhivhr) - jhiuhr) - (/, Vh - Uh)L^{n)2 

= a{uh,Vh - Wh) + b[vh - Wh,Ph) -~ {f,Vh - Wh)L2(ny 

+ a{uh,Wh - Uh) + ihiwhr) - 3h{uhr) ^ {f,Wh - Uh) mn}^ 

>o. 

Hence {uh,Ph) is a solution of VI^. 

o 

(b) Let (uhjPh) G Vnh X Qh be a solution of Ylh- Taking Uh ± Vh as a test 



function in (4.1), with an arbitrary Vh G Vh, we have 



a{uh,Vh) + b{vh,Ph) = if,Vh)L2(n)2 (Vw^ G Vh). 



(4.10) 
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Therefore, since {vh G Vnh \ (vhr, 'nh)Ah = (Vr/;i e Ah)} = Vh, the inf-sup condition 
given in Lemma |3.2| asserts the unique existence of Xh € A;i such that 

a{uh,Vh) + b{vh,Ph) + ivhT,Xh)A,, = (/,w/i)L2(n)2 (Vw/i G Vnh)- 



Combining (4.11) with ( |4.1| , we obtain 

{VhT - Uhr, A/Oa,. < jhivhr) - jhiuhr) 

which gives, by a triangle inequahty, that 

{Vhr - Uhr, Xh)Ah < jhiVhr " Uhr) 



(Vw/i &Vnh), 
{VVh € Vnh)- 



From (4.13) together with Lemma |3.3[ i), we deduce 

iVh, Xh)Au < jhiVh) {yvh G A/i). 

Hence Lemma 



(4.11) 
(4.12) 
(4.13) 
(4.14) 



3.4 



iii) imphes that Xh G Ah, and (4.4) is estabhshed. It remains 
only to prove (4.6). Taking Vh = in (4.12), we have jhiuhr) < {uhT,Xh)Ah- This 
implies (4.6) by Lemma 3.4 ii). Therefore, (uh,Ph, Xh) is a solution of Problem 

o ^ 

(c) Let {uh,Ph, Xh) G Vnh X Qh X A/i be a solution of Problem VEh- Then we 
see that Uh G Vnh,a from (4.5), and that 

{uhT, Xh)Ah 



jhiuhr) 



(4.15) 



from (4.6) combined with Lemma [3.4[ ii). It follows from (4.4) and (4.15) that 

a{Uh,Vh - Uh) + jh{vhr) - jhiUhr) " (/, Vh - U/i)l2(o)2 
= -b{vh ~Uh,Ph) - (vhr -Uhr,Xh)Ah + jhiuhr) - jhiuhr) 
= - (vhr - UhT, Xh)Ah + jhivhr) - 3h{uhr) 
=ih{Vh) - {vhT,Xh)Ah 
>0 {^Vh G Vnh,a), 



where the last inequality holds by Lemma 3.4[i). Hence Uh is a solution of Problem 
VI/i.o-. This completes the proof of Theorem 



4.1 



□ 



4.2. Error analysis. Before presenting the rate-of-convergence results, we state 
the following: 

Proposition 4.1. Let {u,p) be the solution of Problem VI and {uh,Ph), that of 
Problem Ylh for h > 0. Then, 

(i) it holds that 

\\uh\\m{ny < ll/llL2(o)2/a. (4.16) 

o 

(ii) for every Vh G Vnh o,ud qh G Qh, it holds that 

a\\u - < aiu -Uh,u- Vh) + b{uh - u,p- qh) + b{vh - u,ph - p) 

+ {o-T,VhT -Ur)L2(ri) +i{uhT) - jhiuhr) + jhivhr) - i{Ur). (4.17) 

o 

(iii) for every qh G Qh, it holds that 

\\p~Ph\\L^(n)^ < (^1 + ^^ ||p-g/.|lL2(o)2 + i^||u-u,,|lffi(o)2. (4.18) 



Proof, (i) Since Uh is the solution of Problem Ylh,c! by Theorem 4.1 4i), it satisfies 



(4.7). Hence Korn's inequality (2.10), together with the positiveness oi jh, gives 

"II'"/i|Ihi(o)2 < a{uh,Uh) = {f,Uh)L^{ny - jhiuhr) 

< {f,Uh)L^(<n)^ 

< \\f\\L2{ny\\uh\\L^(ny, 
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which imphes (4.16). 

(ii) Let Vh G Vnh and qh G Qh be arbitrary. We begin with the following equality: 

a{u -Uh,u- Uh) = a{u - Uh, u-Vh)~ a{u, Uh - u) - a{uh,Vh - Uh) + a{u, Vh-u). 



We bound from above the second term of the right-hand side by ( 2.16 1 with v = Uh, 
the third one by (4.1) with Vh itself, and rewrite the fourth one by (2.12) with 
V — Vh — u. Consequently, 

a{u -Uh,u- Uh) 

<a{u ~Uh,U~ Vh) + b{uh - U,p) + jiuhr) - jiUr) - if, Uh - u)L2(n)2 

+ b{vh - Uh,Ph) + ihivhr) - jhiuhr) ^ (/, Vh " U/i)L2(f2)2 
- b{vh ~U,p) + {(Tr,VhT - Ur)L^ri) + if: - '")L2(a)2 

=aiu - Uh,u- Vh) + b{uh - u,p - Qh) + bivh - u,ph - p) 

+ icrT,Vhr - Wr)L2(ri) + jiuhr) - jhiuhr) + jhivhr) - 

Combining this with Korn's inequality ( 2.10), we conclude ( 4.17[ ). 

(iii) Taking u ± u as a test function in (2.161, with an arbitrary v G iJQ(ri)^, 
gives 

aiu,v) + hiv,p) = {f,v)L-(n)- (Vi; G HH^)^). 
On the other hand we know that (4.10) holds, and therefore, by subtraction we 
obtain 



Ph 







i^Vh G 



aiu-Uh,Vh) + b{vh,p- 

o 

Now let qh G Qh- It is clear that 

Wp-PhWh'^in) < lb - g/i||L2(n) + hh-PhWh-^in)- 
By Lemma 3.1^1) together with (4.19), we have 



(4.19) 
(4.20) 



-Ph\\ < sup 

Vh&Vh 



bivh,qh~Ph) 



sup 



bivh,qh -p) +bivh,p-ph) 



sup 

vheVh 



\\Vh\\H^n)^ 
bivh, qh - p) - aiu - Uh, Vh) 



\Vh\\H^fl)^ 



lb/i||_ffi(0)2 

< ll^ll Wp-QhWL^n) + l|a|| \\u ~ '"^11^1(0)2 



The desired inequality (4.18) follows from (4.20) and (4.21). 



(4.21) 

□ 



We are now in a position to state the primary result of our error estimates, 
assuming only the regularity of the exact solution. 

Theorem 4.2. Let (u,p) be the solution of Problem VI and iuh,Ph) be that of 
Problem Vlh for < h < 1. Suppose g G C\Ti) and {u,p) E H^+'in)^ x H'{n) 
with < e < 2. Then we have 



\\u- Uh\\Hi(ny^ + \\p-Ph\\L^n) < Ch" 



(4.22) 



Proof. We recall the interpolation error estimates (3.1|-(3.3|. 
Taking Vh = X^m, qh = TihP in (4.17) and (4.18), we find that 

"lb - ■"/i|lffi(n)2 < a(u - Uh,u-Ihu) + biuh - u,p - Ilhp) + bilhu ~ u,ph - p) 

+ iar,Ur - iThu)r)L^(ri) + Uiuhr) " jh('"/ir)| 

+ \jhiilhu)r)-j{ilhu)r)\ + \3iiIhu)r)-jiUr)\, (4.23) 

and that 

\\p-ph\\mn) < c(lb-n/ip|ii2(f^) + lb - u,i||^i(Q)2) 

<Cih' + \\u-Uh\\miny)- (4.24) 
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Each term of the right-hand side in (4.23) is estimated as foUows: 
1. 

\a{u- Uh,u-Ihu)\ < \\a\\ \\u - Uh\\HHny\\u - Ihu\\Hi(n)2 < Ch^Wu - Uh\\Hi{ny. 
2. 

\b{ui-,~ u,p-Uhp)\ < \\b\\ \\u- Uh\\Hi{nr\\p-'^hP\\L^n) < C^'ll^^ - "/ilki(a)2- 
3. From ( |4.24| , 

\b{lhu-u,ph-p)\ < \\b\\ \\lhu - u\\H^n)^\\Ph - p\\mn) 
< C(h^' + h'\\u~Uh\\HHny)- 

4. 

{ar,Ur ~ {Ihu)r)^2(^Y'i) ^ I | (r^ ) 1 1 Ur - (Xft u)^ 1 1 ^2 (r^) < C/l3+'. 

5. By Lemma [3^ ii) together with Proposition |4.1[ i), 

iJiUhr) - jh{uhr)\ < Ch^WuhrWjjlf^^^^ < C/l^ 1 1 | l < Ch^. 

6. Since \\Ihu\\Hi{n)^ <\\^hU ~ u\\Hi{n)2 + ||w||ffi(o)2 < C, Lemma [S^ejii) imphes 
7. 

\j{{Ihu)r) ~ j{ur)\ < \\g\\L^(ri)\\ur - {Ihu)r\\mri) < Ch^^" . 
Combining these seven estimates with (4.231, we deduce that 

II" - "'iIIhi(o)2 < C {h^Wu - Uh\\m{n)^ + + h^+'' + h^^ 



- jh{{Thu)r)\ < Ch^\\{Ihu)r\\^i^^^^ < Ch^\\Ihu\\Hi(n)^ < Ch^ . 



Therefore, 



I"" - Uh\\H^{ny < C/i" 



We conclude (4.22) from (4.24) and (4.25) and this completes the proof. 



(4.25) 

□ 



The previous theorem reveals that the rate of convergence is 0(/ii/4) at best even 
when the solution is sufhciently smooth. However, it can be improved if additional 
conditions about the signs of Uf„- and (Ihu)r on Fi are available. To formulate the 
result, we make the following assumptions (Recall Definition 3.1 and see Remark 



4.1): 



(51) (I/iu)r has a constant sign on every side. 

(52) Uhr has a constant sign on every side. 

(53) sgn(u^) = sgn((I,ju)^) on Fi. 



Theorem 4.3. In addition to the hypotheses in Theorem 4.2, we assume g € 
C^(Fi) and that (SI) -(S3) are satisfied. Then we have 

\\u-Uh\\H^n)- + \\p^Ph\\mn)<Ch^'-^'''K (4.26) 

Moreover, if g is a polynomial function of degree <1, we have 

\\u- UhWH^iny + \\p-Ph\\L^n) < Ch\ (4.27) 

Proof. We first verify that (S3) implies 

arii:hu)r + g\ilhu)r\ = a.c. ouFi. (4.28) 

In fact, for each side e € <o';i|ri, if Ur vanishes on a subset of e containing more 
than three points, then the quadratic polynomial {Ihu)r vanishes on the whole e. 
Otherwise, we have \ur\ > a.e. on e; hence we deduce from (2.15), namely, 

a^^Ur + g\ur\ — a.e. on Fi, (4.29) 
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that (It = —gsgn{ur) ~ — .gsgn((I^M)T-) a.e. on e. In both cases, it fohows that 
ariIhu)T + g\{'Ihu)T\ — a.e. on e. Thus (4.28) is vahd. 
It follows from ( |4.28p and ( |4.29[ ) that 

Therefore, taking Vh = 1-hU and qh = H-hP in (4.17) gives 
a\\u - < a{u - Uh,u-Ihu) + b{uh - u,p - Xlhp) + b{IhU - u,ph ~ p) 

+ \i{uhr) - ih{uhr)\ + \ih{{Ihu)r) - ]{{IhU)r)\- (4.30) 

Let us give estimates for each term on the right-hand side. We can evaluate the 
first three terms by the same way as in the proof of Theorem |4.2| By Lemma |3.7[ 
the fourth and fifth terms are estimated as 

\j{uhT) - jh{uhT)\ < Ch^WuhrWh^iVi) < Ch^\\uh\\Hi(ny < Ch^ , 

\jh{.{Xhu)r) - ]{{IhU)r)\<Ce\\{IhU)r\\L-[r^) < CK^WIhuWh^U)- < Ch\ 

Consequently, we obtain 

II" - UhWHi^ny <C{h'\\u- UhWn^ny + h^' + h^) , (4.31) 

which leads to 

The estimate for \\p — Ph\\L2{n) is similar to the proof of Theorem 4.2 and then, 

( |426l ) follows. 

Finally, if g is affine then the fourth and fifth terms in (4.30) vanish exactly, 
according to Lemma |3.7[ Hence we have 



II" - "'i|Ihi(o)2 < C {h^Wu - Uh\\H\n)^ + h'^") 



instead of (4.31), from which (4.27) follows. 



□ 



Remark 4.1. Conditions (S1)--(S3) are not so artificial. Assume that u, the velocity 
part of the solution, is continuous on Q and that the isolated zeros of on Fi are 
contained in Ti_h- If we make h sufficiently small, then we see that (SI) and (S3) 
are satisfied. Therefore, since Theorem 4.2 implies Uhr Ut in i?5(Fi), we can 



expect (S2) to also be valid; however, its rigorous proof is not easy. 

4.3. Numerical realization. We propose the following Uzawa-type method to 
compute the solution of Problem VE/j (therefore. Problem VI/j) numerically. 

Algorithm 4.1. Choose an arbitrary X^^'' e Ah and p > 0. Iterate the following 
two steps for A; = 1, 2, • • • : 

Step 1 With X^i^^ known, determine (u;f\p|f^) € x Qh by 



Step 2 Renew A|f+^^ G A,, by 



A 



(fc+i) 



Proj^JAi'=)+K;0 



(yvheVnh), (4.32) 
(yqheQh). (4.33) 

(4.34) 



Remark 4.2. (i) The unique existence of iu\^\p'h^) satisfying (4.32) and (4.33) is 
guaranteed by the inf-sup condition mentioned in Remark |3.1[ 



(ii) We can regard (4.34) as an approximation of 

Xh = PtoJa^ (A/i + pUhr), 



(4.35) 



which is equivalent to (4.6) by Lemma 3.4111 
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Theorem 4.4. Let (uh,Ph, ^h) be the solution of Problem VE^. Under the same 
notation as Algorithm |4.H there exists a constant po > independent of h such 
that if p satisfies < p < po, then the iterative solution {u')^ \ p\^ \ X'^j^^ ) converges 
to {uh,Ph, Xh) in H^{^})^ X L^{fl) x Ah, as k oo. 



Proof. Subtracting (4.32) from (4.4 1 with test functions in Vnh,a, we obtain 



(4.36) 



Uh £ Vnh,a and apply Korn's inequality (2.101 to 



In particular, we take 
obtain 

(^hr -W'ir,Ajf^-A,OA^ ^ ~a{u[^'' -Uh,u^^^ -Uh) < -a||u*f^-u/i||^i(Qp. (4.37) 



(4.38) 



Next, we note that Proj^^ given in (3.6 1 satisfies 

as a result of a general property of a projec tion op erato r. It f ollow s from (4.38) 
with ph = ^ + pu\^^ and r]h = Xh + puhr-, (|4.34|, (|4.35[), and (IZstI that 



l|A 



(fe+i) 



Xh\\l,<\\Xi^^-Xh+p{u):>~Uhr)Wt 



(fc) 



= I|A 
< IIA 



(fe) 

h 

(fc) 



A/.IIL + Mu^hr - ^hr, Xr - A,)a, + p^u'C; - UHrWl, 



(fc) 



(fe) 



Therefore, since \\u''^^ - "ftrlU,, < C'll^J^V - Uhr\\L2(ri) < C\\u''^' - Uh\\Hi(ny m 
view of Lemmas 3.6 |i) and 2.4 'i), we obtain 



,(fe) 



llAi^^) - A.lli < ||aW - A.lli - (2«p - Cp')\\ui'^ - u.ll^.(^). 
and thus 

(2ap - C7p2)hf - UhWl.^^y < l|A,?^ - A.lli,^ - |lA|f+^' - A,||i, 



(4.39) 



On the other hand, by virtue of Lemma 3.3 ^i)(ii), we can choose Wh £ V^^a such 
that WhT = A|f "* — Xh on Fi and 

Ak 



\\wh\\m(nr < C\\Xr - Xh\\m/Hv,) < Cih)\\x['^ - A„||a,, 

where the constant C{h) concerns the equivalence of the norms on the finite dimen- 
sional space Ah- Hence, it follows from (4.36) with Vh = Wh that 

ll^i''^ - ^hWl^ = {whT,Xh ^ - Xh)A^ = -a(Ai''^ - Xh,Wh) 

< l|a|l W^h ^ - Uh\\H^n)A\wh\\H^ny 



<C{h)\\u['' -UhWHHnyUr - XhU, 



(fe) 



so that 



Cih)\\ui^' - UhWminr < -||Ar ->^h\\A 



(fe) 



(4.41) 



Since the constant C in (4.39) is inde pend ent of p (a nd even of /i), if we choose 
< p < Po := ^ then it follows from |4.39| and ( |4.4l| that 



l|A,i -A/i||a, 



< 



/ 2ap-Cp^ (k) . n 



where we may assume C{h)^ > 2ap ~ Cp^ (if not, take C{h)^ — o? jC). Conse- 
quently, we conclude 



IIA 



(fe) 
h 



X 



/i A, 



< 1 



lap - Gp^ 
C(hY 



Wx'^h^ -Xh\W^^ (fc^oo). 
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Then, from (4.40) it also follows that u-y^^ Uh in H^{VtY as fc — > oo 



Finally, subtracting (4.32 1 from (4.4) with test functions in Vt gives 



P\\Ph' -PhWh^iO.) < sup 



a{uh - ul^\vh) + b{vh,Ph - P^h'') = {yvh eV/i). 

a(u]^ - Uh,Vh) 



Therefore, by Lemma 3.1 4) we have 

u (fc) 



Ph 



\Vh\\H^(ny 



< Ikll \\Uh^ - "/illifi(n)2 
^0 (fc ^ oo). 



= sup 



\vh\\H^(ny 



This completes the proof. 



□ 



5. Discretization of the Stokes problem with LBCF 

5.1. Existence and uniqueness results. Approximate problems to Problem VI 
(therefore. Problem PDE) in the case of LBCF are as follows. 

Problem 5.1 (VI^). Find {ufi,ph) € Vrh x Qh such that 

a{Uh,Vh-Uh)+b{vh-Uh,Ph)+jh{Vhn)-jh{Uhn)>{f,Vh-Uh)mn}^ (yvhEVrh), (5.1) 

b{uh,qh)^0 (yqhEQh)- (5.2) 

Problem 5.2 (VI^.u)- Find Uh E Vrh,a such that 

a{uh,Vh - Uh) +jh{vhn) - jh{uhn) > {f,Vh - Uh)L^(ny (Vu/,. e K^.^)- (5.3) 
Problem 5.3 (VE^). Find {uh,Ph,^h) G Vrh x Qh x A/,, such that 

aiuh, Vh) + b{vh,Ph) + {vhn, A?i)Ah = (/, w/i)L2(f2)2 {yvh e Vrh), 
b{uh,qh)^0 (yqh^Qh), 
(u/m, M/i - A/i)a^ < {Vfih e A,i)- 



(5.4) 
(5.5) 
(5.6) 

Further- 



Theorem 5.1. (i) Problem Vlh.a admits a unique solution Uh G V^t/i, 
more, it satisfies the following equation: 

a{uh,Uh) + jh{uhn) = {f,Uh)L^ny- (5-7) 

(ii) Problems Ylh and YEh are equivalent in the following sense. 

(a) // {uh,Ph) S Vrh X Qh is a solution of Problem VI^, then there exists a unique 
Xh G A/j such that {uh,Ph, Xh) solves Problem VE^. 

(b) If {uh,Ph, Xh) G Vrh X Qh X Ah is a solution of Problem VE/i, then {uh,Ph) 
solves Problem VI/j. 



Proof. Statements (i) and (ii)(a) are proved by the same way as Theorems 4.14) 



and (ii)(b), respectively. We obtain (ii)(b) by a discussion similar to that in the 
proof of Theorem |4.1| [^ii) (c) . □ 

Remark 5.1. It is clear that if {uh,Ph) is a solution of Problem VI/j, then Uh is 
that of Problem Vlh,cr- However, the "converse" is no longer true. In fact, the 
pressure, especially its additive constant, need not to be uniquely determined even 
if the solution Uh of Problem Ylh,a is given. This situation is quite different from 
the case of SBCF. 



The next theorem guarantees the existence of a solution of Problem Ylh- 
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Theorem 5.2. (i) There exists at least one solution of Problem VE/j, and the 
velocity part is unique. 

(ii) If {uh,PhT^h) ciiT'd {''J-inPhT ^h) '^^^ ^'""^ solutions of Problem WEj^, then there 
exists a unique G R such that 

Ph=Pl + 6h and K{M)^\l{M) + ^ {^MeVi.h)- (5.8) 

(iii) Under the assumptions in (ii), if we suppose Uhn 7^ onVi, then 5h — 0. 
Namely, a solution of Problem VE/i is unique. 



Proof, (i) The uniqueness of the velocity is obvious; see Remark 5.1 



To prove the existence, let Uh £ Vh.a be the solution of Yl^.a- Taking Uh ± Vh 



as a test function in (5.31, with an arbitrary Vh €Vh H Vrh.a, we obtain 

o 

a{Uh,Vh) = (/,W/i)L2(f2)2 (VWft €Vh n Vrh,a)- 

o 

Hence, from Lemma 3.1 'i), we deduce the unique existence of G Qt such that 



a{uh,Vh) + b{vh,pl) = (/, Wh)L2(a)2 {\/vh €14). 

o 

Moreover, noting that {vh € Vrh \ ivhn,ilh)Ai, = (Vyy^ G Ah)} =Vh, it follows from 



Lemma 



3.3 



iii) that there exists a unique A^^ G A/j satisfying 



a{uh,Vh) + b{vh,pl) + {vhn,^h)Ah = {f,Vh)mny C^Vh e Vrh)- (5.9) 



For every rjh £ Ah Ci LQ(ri), by Lemma 3.3 'ii) we can choose Wh G Vrh.a- such 
that 

Whn=Vh on Fi. 



Hence, taking Vh = Wh in (5.31 and ( |5.9[ ), we obtain 

a{uh,Wh - Uh) +jh{Vh) - jh{uh„) > if,Wh - Uh)L^(nY 
a{uh,Wh) + [r^hAD^h = {f^Wh)L-^(ny. 



(5.10) 
(5.11) 



Substituting (5.7) and (5.11) into (5.10) gives 



(%, A^)a, < JhiVh) {yVh (^AhH LliT,)). 

We apply Hahn-Banach's theorem to deduce the existence of some A/j G Ah such 
that 



{■nhAh)Kh < jhivh) (y-nh G Ah). 



Therefore, Lemma 



3.4 



ii) implies that \h G Ah. Furthermore, since A^ satisfies 
iVh, Xh - A^)a, = {yr^h eAhH LliVi)), 
it follows from Lemma [3^ iv) that there exists some 6h G R such that 

Sh 



\h{M) = \1{M) + 



9{M) 



(VMGFi,^) 
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Thus, from Simpson's formula and (5.9), we obtain 



--a{uh,Vh) +b{vh,Ph.) + («/m,A")A^ 



6 

i—1 



(vhnA 



4v 



hn,i-\ 



^a{uh,Vh) + b{vh,Ph + ^h) + {vhn,^h)A„ (Vw?i e Vrh)- 



(5.12) 

(5.13) 
(5.14) 
(5.15) 



This estabhshes (5.4) if we define 



pI 



5h. 



Equation (5.5) obviously holds because £ Vh^cr- It remains to show (5.6), 
which is equivalent to {uhm^h)^^ — ih{uhn) by Lemma 3.4 li). This is indeed 
obtained from (5.15) with Vh = uu and (5.7). 

(ii) Let {uh,Ph,^h) and {uh,p*i^,X*i^) be two solutions of Problem VE;i. Because 
the uniqueness of the pressure up to additive constants is shown in the proof of (i), 
there exists a unique constant 5^ such that pf^ — p'^ + S^- 

Since {uh,Ph, ^h) and (u/i,p^,A^) satisfy (5.4), subtracting the two equations 
and calculating in a manner similar to (5.12 )- ^5.15 ), we obtain 

{Vhn/^h)Ah (VWh e Vhr)- (5.16) 

Here, A,. G A,, is defined by Xh{M) = Xh{M)- XUM) - 5h/g{M) for M £ Tim- It 



follows from ( 5.16 ) toge ther with Lemma 3.3 4) that (77/;, A/Ja^ = for all rjh G A/j 
Hence Xh — 0, and (5.8) is proved. 

(iii) The assumption Uhn 7^ implies that either of the following is true: 

o 

(a) There exists A/GFi^/i such that Uhn{M) > 0. 

(b) There exists MGLi./i such that u/i„(M) < 0. 
Since Uf^ G Vrh,a^ we see that 

Uhn ds — div Ufi dx = 0. 



Therefo re, b oth (a) and (b) above must hold true. Then, from (5.6) combined with 

o 

Lemma 3.4 4i), we find some M,N G Ti^h satisfying Xh{M) — 1 and Xh{N) — —1. 
Consequently, the additive constant Sh appearing in (ii) cannot att ain any value 
except because Xh E Ah- This completes the proof of Theorem 



5.2 



□ 



5.2. Error analysis. Let us begin with the following analogue of Proposition |4.1[ 

Proposition 5.1. Let {u,p) be a solution of Problem VI and {uh,Ph) be that of 
Problem VI/,. for h > 0. Then, 

(i) it holds that 

< ||/||L2(n)2/a. 

(ii) for every Vh G Vrh o,nd qh G Qh, it holds that 

a\\u - Uh\\H^Q^2 < a{u -Uh,u- Vh) + b{uh - u,p- qh) + b(vh - u,ph - p) 

+ {(^n.Vhn - W„)i2(ri) +j{uh„) ~ jh{uhn) + jh{vhn) ~ j (Un) ■ (5.17) 
o 

(iii) for every qh G Qh, it holds that 



\\P -Ph 



7JU^(o)< (^i+^^y\p'-qh\\LHn)-+^^\\u-Uh\\HHny, (5.18) 
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where p° ^ p - {p, I) ^2 (^^^^ / \n\ and pi = ph - {ph, l)L^n)/\^\- 
(iv) for every e Qh, it holds that 

Wp-PhWmn) < (l + ^) \\p-<lh\\mn) + ^^\\u-Uh\\HHnr-+C. (5.19) 

Proof. Statements other than (iv) can be proved by the same way as Proposition 



4 J] To show (iv), we let Qh G Qh- It is clear that \\p-Ph\\L^(n) < \\p - <lh\\ l-^ (n) + 
Qh ^ Ph\\L'^{n)- To bound the latter term, we deduce from Lemma 3.1 together 
with ( |2.18| and ([5^, that 

All II ^ b{vh,qh-Ph) b(vh,qh-p) + b{'"h,P-Ph) 
P\\qh-Ph\\L^{n) < sup — — = sup n — n 



sup 



b{vh,qh -P)- a(" - Uh,Vh) + {(Tn,Vhn)L-^(Tt) + [Vhn, K)Ah 



||«/l||_ffl(fl)2 

< \\p - qh\\L^(n) + \\a\\ \\u - ?i/i||ffi(o)2 + C. 



Here, to derive (5.20), we have used the estimates 



(5.20) 

< C\\vh\\H\n)^, 



which are obtained from Lemmas 2.3 ^i), pO| i), and 3.6 i). The desired inequality 
(5.19) immediately follows from (5.20). □ 



Remark 5.2. lip — ph e io(^) Ih & Qh, we can ta ke C = in (5.19) according 
to the equality p — ph = p'^ ^ Ph combined with (5.18). 

We state the rate of convergence result for the case of LBCF, which is not better 
than that of SBCF because of the influence of an additive constant of the pressure. 

Theorem 5.3. Let {u,p) be a solution of Problem VI and {uh,Ph) be that of Prob- 
lem VI/i for < h<l, and suppose {u,p) G H^+''{n)^ x iJ'(f2) with < e < 2. 
Then, 

||"-^h||ifi(o)2 + ||/-p?JU2(o) <C/i-'"{4'i>, (5.21) 
where p° =p - (p, l)i2(n)/|f^| and p^ = pu - {phA)L'^(n)IM- 

IhU and qii 



Proof. Let us take Vh 



H-hP in (5.17) and bound from above each 



term on the right-hand side. By (5.19), we have 

\b{lhu - u,ph -p)\ < \\b\\ \\lhu - u\\Hi(n)^\\Ph - p\\L^-{n) (5.22) 
< Cih'Wu - UhWminv + + ^')- (5-23) 

For the other terms, we employ the same estimates as those in the proof of Theorem 
EH Then it follows that 



which implies 



\\p"-pl\\LHn) 



W/!.||//l(f2)2 



< Ch" 



(5.24) 



(5.25) 



< C'/i™'"^2'4} follows from (5.18) and (5.25) and this completes the 



proof. 



□ 



Remark 5.3. (i) If we assume, in addition, that p — Ph & -^o(^) then we can 
establish the result of 0(/i™"'t'^'j'^). Moreover, as in the case of SBCF, under 
suitable conditions regarding the signs of Uhn and {Ihu)n on Fi, it can be improved 
to 0(/i"i'"{''i>), or even 0(/i') if g is affine. 
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(ii) When the uniqueness of Problem VI holds, we can obtain a strong conver- 
gence result for the error of the pressure including the additive constant. In fact, 
the uniform boundedness of ph in L^(fl) gives a weak convergence limit p for some 
subsequence p/j/. Since u/j — u in H^(fl)^, we have 

\jh{uhn) - < \jh{uhn) - j(u/m)| + \jiuhn) -j{Un)\ ^0 {h ^ 0). 



Therefore, taking the limit /i' — > in ( |5.1| -( 5.2 ), we find that {u,p) is a solution of 
Problem VI. Hence p = p, and from ph' — ^ p and pi — > we conclude the strong 
convergence of the whole sequence. 

5.3. Numerical realization. Based on Problem VE/i, we propose the following 
Uzawa-type method to compute the approximate solution numerically. 

Algorithm 5.1. Choose an arbitrary aJ^^-* G A/j and p > 0. Iterate the following 
two steps for /c = 1, 2, • • • : 

Step 1 With x\^^ known, determine iu\^ \pl^'') G Vrh x Qh by 



a{u[^\vh) + b{vh,Ph^) = {f-,Vh)v^(ny - {vhnAh^)Ah 



step 2 Renew A^f € A,, by 



Ak+l) 



= Pro,iA.(Af 



i^vueVrh), (5.26) 
(Vg/,eg„). (5.27) 

(5.28) 



(k) (k)\ 

Remark 5.4. The unique existence of [uj^ ,p]i ) satisfying (5.26) and (5.27) is 



guaranteed by the inf-sup condition in Lemma |3.2[ 

Theorem 5.4. Under the same notation as Algorithm |5.1[ there exists a constant 
Po > independent of h such that if p satisfies < p < pq, then {u^^\ p''y^\ \^^^ ) 
converges to some solution of Problem VE^ in i?^(J7)^ x L^(r2) x A^, as k oo. 

Proof. First wc show the boundedness of the sequence {{u\^\p\'^\ '^\!^^)}k- In fact, 
taking Vh = in (5.261, we find from Korn's inequality (2.10) that 



II (fe)l|2 



which gives ||ffi(si)2 < C Here, to derive the last line, we have used 



<ll/llL^(n)HI"i'^IU^(o)^ + l("i';\Ai'=^)Aj 



(fe) ,(fe)^ 



|(uit\Af )aJ <j,(<';0 <C|l<„^IU^(ro <C^II<^1ifno)^ 



,ik) 



(fc)| 



(5.29) 

which is obtained from Lemmas |3.4[ i), |3.6| [i), and |2.3| ^i). Then Lemma 3.2 together 
with (5.261, implies that 

All (fe)ll / Hvh.Ph^) 

^m'\\LHny< sup 1^ — ^ — 



sup 



l|w/i||ffi(n)2 



< 11/1^2(0)2 + ||a|l \\u\^^\\m{n) 

< c, 



sup 



where {{vim, '^|f^)Afc| hi the third line is estimated in a manner similar to (5.29). It 



is clear that {X\^^}k is bounded because X\^' G Ah- 



(k) 
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Therefore, we can extract a subsequence {{u\^ \Ph ^ '')}k' converging to some 
element {uh,Ph,^h) G Vrh 'xQh'x A^. Making k — k' and fc' — > cx) in (5.26)-(5.28), 
we obtain 

a{uh,Vh) + b{vh,Ph) + {vhnAh)\u = if^Vh)L^(n)'^ (Vw^ G Vrh), 
b(uh,qh) = [yqh^Qh), 
Xh = Proj^^ {Xh + puhn)- 



Consequently, since the last equation is equivalent to (5.6) by virtue of Lemma 



3.4 li), we see that {uh,Ph, Xh) is a solution of Problem VE/j. 

It remains only to prove that the whole sequence converges to (uh,Ph, Xh)- Sub- 
tracting (5.26) from (5.4), we have 

a{uh- u\^\vh) + b{vh,Ph - Ph^) + [vhnAh - >^h^)AH ^ (Vwft e K/i)- (5.30) 



(k) 

In particular, if we take Vh = Uh — uj^ , then 

Therefore, it follows from a general property of the projection operator Proj^^ that 

||A.-A('=+i)||i<||A,-Ai'=)+pK„-.W)||2^ 



, (fc) (fc)x 

-a(uh ~ u\' ,Uh- ul')- 



l|Ah - Ai'^'lli,^ + {Cp' - 2ap)\\uh - ui'^\\l,^,,y, (5.31) 



,W||2 



< 



where we have applied Korn's inequality (2.10) and the estimate — u[^J^\ 



< 



C\\uhn — u'"j^n\\L^Cri) — C^\\''^h — u\i'\\H^{ny to derive the last line. Since the constant 
C in (5.31) is independent of p (and even of h), we choose < p < po := ^ to 
obtain 



,(k) 



0<\\uh-uf^\\l,^,,^. 



o < 



1 



2ap - Cp^ 



l|A.-Af ||L-||A.-Ar^^| 



(fc+l)||2 



A, 



(5.32) 



Hence the sequence {\\Xh — X^l^\\\^}k is decreasing. Noting that a decreasing se- 
quence in R, bounded from below, converges to its infimum and that inf^ \\Xh — 
AJi llAfc = by construction, we conclude that X\ ^ Xh in A?i as fc -> oo. From 
this and (5.32), it also follows that u'^'' — > Uh in H^iVlf'. Finally, from the inf-sup 
condition given in Lemma 3.2 combined with (5.30), we have 

All {k)u ^ b{'"h,Ph-Ph^) 

PWPh-pi'h'^n) < sup — ^— ^ 

.ik] 



sup 



'a{uh 



,Vh) - {vhn,Xh - A|j'''')a^ 



< \\u^h^ - Uh\\m(i^)^ + ||Ai'°^ - A^IIa;. sup 

< \\a\\ \\u[^^ - UhU^n)- + C\\X[^^ - Xhh, ^0 (fc ^ oo). 
This completes the proof. □ 



Remark 5.5. (i) The resulting solution of Problem VE/i as the limit of (u\^\ p\^\ X^''^ 



h )■• 



especially its additive constant of the pressure, may depend on a choice of < p < 
Pq or that of the starting value X^^\ However, if Uhn ^ on Fi, and hence the 
uniqueness of the solution of Problem VE^, is valid, then it is obviously independent 
of them. 
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(ii) Contrary to the case of SBCF, it is difRcult to prove an exponential con- 
vergence of the iterative solution (u'^^ , p^^^ , a|/^' ) because we do not know whether 



A/i — \^^^ e iQ(ri), which is necessary to deduce an extension of \h — X)^' to Vrh,(7- 

6. Numerical examples 

We assume fl — (0, 1)^, the boundary of which consists of two portions Tq and 
Fi given by 

To = {(0,y) I < 2/ < 1} U {{x,0) | < x < 1} U {{l,y) | < y < 1}, (6.1) 
Ti = {(x,!) |0 < X < 1}. (6.2) 

In particular, the set of extreme points is Fq n Fi = {(0, 1), (1, 1)}. For the trian- 
gulation 5^ oi Q, we employ a uniform N x N Friedrichs-Keller type mesh, where 
N denotes the division number of each side of the square fl. 
Let us consider 

= 20x^(1 -x)^yil-y){l~2y), 
= ~20x{l - x){l - 2x)y^{l - y)2, 
= 40a;(l -x)il~ 2x)yil - y){l - 2y) 

+4(6x5 - 15x4 + I0x^){2y - 1) - 2, 

which turns out to be the solution of the Stokes equations under the adhesive 
boundary condition for u = 1 and / given by 

/i(x,y) = 0, 

/2 (x, y) = 120(2x- 1)?/2 (1 - y)2 + 80x(l -x) (1 - 2x) (6^^ -6y + 1) +8(6x5 _ ^^^i ^ -^Qx^ . 
By direct computation, we have 



(fe) 



'ui{x,y) 
U2{x,y) 
P{x,y) 



(6.3) 



rnax |cr^| 

Ti 



max 20x^(1 

0<a;<l ' 



xf 



4 - 1-25' 



max|cr,J= max — 4(6x 

Pi 0<2;<1 ' 



15x'' + 10x^) + 2| =2. 



(6.4) 
(6.5) 



Now, if we impose SBCF or LBCF on Fi, with g being constant, instead of the 
adhesive boundary condition, then in the case of SBCF, we find that 




(6.3 1 remains a solution. 



(6.3 1 is no longer a solution and a non-trivial slip occurs. 



and in the case of LBCF, 



5>2 
g<2 



(6.31 remains a solution. 



(6.3) is no longer a solution and a non-trivial leak occurs. 



We indeed observe some of the abovementioncd phenomena in our numerical com- 
putation, as indicated in the plots of the velocity field shown below in Figures 6.1 



and 6.2 In addition, we find that the bigger (resp. smaller) the threshold g of a 
tangential or normal stress becomes, the more difRcult (resp. easier) it becomes for 
a non-trivial slip or leak to occur, which is in agreement with our natural intuition. 



Next, we consider the behavior of th e Lagrange multiplier A^. It follows from 



(4.6) or (5.6) together with Lemma 

{\K{M)\<1 
\Xh{M) 



3.4 



1 or 



ii) that for each M S Vi^h 

if u^(M) =0, 
if uh{M)^Q, 



1 



(6.6) 



which is observed by comparing the result of Table [6TT] with Figure 6.1 or 6.2 



In Table 6.1 we also see that if any leak does not occur, then the choice of the 
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(a) g = 0.1 



(b) 5 = 0.8 



(c) g = 2.0 



Figure 6.1. Solution velocity field of the Stokes equations with SBCF 





(a) g = 0.1 



(b) 5 = 1.2 



(c) g = 3.0 



Figure 6.2. Solution velocity field of the Stokes equations with LBCF 
Table 6.1. Values of the Lagrange multiplier A/i on Fi 





SBCF 


LBCF 


9 


0.1 


0.8 


2.0 


0.1 


1.2 


3.0 


3.0 


P 


1000.0 


50.0 


3.0 


20.0 


30.0 


2.0 


2.0 




0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.2 


X 








0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.1 


-1.0 


-0.26 


-0.09 


-1.0 


-1.0 


-0.63 


-0.43 


0.2 


-1.0 


-0.90 


-0.25 


-1.0 


-1.0 


-0.57 


-0.37 


0.3 


-1.0 


-1.0 


-0.42 


-1.0 


-1.0 


-0.45 


-0.25 


0.4 


-1.0 


-1.0 


-0.55 


-1.0 


-0.83 


-0.25 


-0.05 


0.5 


-1.0 


-1.0 


-0.60 


-0.06 


-0.06 


-0.02 


0.18 


0.6 


-1.0 


-1.0 


-0.55 


1.0 


0.67 


0.22 


0.42 


0.7 


-1.0 


-1.0 


-0.43 


1.0 


1.0 


0.43 


0.63 


0.8 


-1.0 


-0.94 


-0.26 


1.0 


1.0 


0.58 


0.78 


0.9 


-1.0 


-0.26 


-0.09 


1.0 


1.0 


0.66 


0.86 


1.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


^itr 


4 


18 


29 


21 


12 


29 


30 
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Table 6.2. Convergence behavior of \\uf, — Ui.ct||/fi(o)2 and \\pii — PidW (n) 



N 


SBCF 


LBCF 




-error 


rate 


L'^-crror 


rate 


iJ^-crror 


rate 


L^-error 


rate 


10 


1.6E-2 




1.6E-2 




1.4E-2 




1.3E-2 




12 


l.lE-2 


1.9 


l.lE-2 


2.0 


l.OE-2 


1.8 


9.7E-3 


1.8 


15 


7.0E-3 


2.1 


6.3E-3 


2.5 


6.4E-3 


2.0 


5.8E-3 


2.2 


20 


3.9E-3 


2.0 


3.5E-3 


2.1 


3.7E-3 


1.9 


3.3E-3 


1.9 


24 


2.6E-3 


2.1 


2.7E-3 


1.3 


2.5E-3 


2.2 


2.2E-3 


2.2 


30 


1.7E-3 


2.0 


1.5E-3 


2.6 


1.6E-3 


2.0 


1.5E-3 


1.9 


40 


9.0E-4 


2.1 


8.5E-4 


2.0 


8.4E-4 


2.2 


8.0E-4 


2.2 



starting value Ajj^'* affects the resulting Uniit (A^^'' = 0.2 in the last column imphes 



(1) 



that \\^\m) = 0.2 for each AI ^Tiji), whereas changing the value of p does not 



cause such phenomena. Here, all the computations shown in Figures 6.1-6.2 and 



Table 6.1 are performed for A'^ = 10 until the stopping criterion 



,(fe) 



< 10- 



(6.7) 



is satisfied in Algorithm 4.1 or 5.1 The number of iterations required to attain 



(6.7) is denoted by /citr. 



Finally, we evaluate the error between approximate solutions and exact ones as 
the division number A'' increases, when g = 0.8 and g = 1.2 for the case of SBCF and 
LBCF, respectively. Since we do not know the explicit exact solutions, we employ 
the approximate solutions with N — 120 as the reference solutions (wrefiProf), and 
numerically calculate \\uh — Mref||ffi(n)2 and \\ph — Pi-ci\\L^(n)- Here, the additive 
constants of phS are chosen such that p/j(0, 0) = pref(0, 0). Then, as Table 6.2 
shows, we can observe the optimal order convergence 0{h^) for both SBCF and 
LBCF. 



7. Conclusion and future works 

A finite element analysis using the P2/P1 element to the Stokes equations under 
SBCF or LBCF is examined. We have proved the existence and uniqueness (partial 
non-uniqueness) results and established the convergence order 0{h^/'^) as error 
estimates for appropriately smooth solutions; sufficient conditions to obtain the 
optimal order 0{h^) are also presented. To compute the approximate solution, 
we have proposed an iterative Uzawa-type algorithm. We have applied it to some 
examples and numerically observed the convergence order of 0{h?). 

In a future study, we would like to extend our theory to a more general situation, 
for example, a smooth domain without corners, nonlinear Navier-Stokes equations, 
a case in which SBCF and LBCF are imposed simultaneously, or a time-dependent 
problem. 
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